To appear in Ap.J. 


The Micro-Arcsecond Scintillation-Induced Variability (MASIV) Survey III. 
Optical Identifications and New Redshifts 

Tapio Pursimo 

Nordic Optical Telescope, Apartado 474, 38700 Santa Cruz de La Palma, Spain 

tpursimoOnot . iac . es 
Roopesh Ojha 

NVI Inc./ U. S. Naval Observatory, 3450 Massachusetts Ave NW, Washington DC, USA 

David L. Jauncey 

CSIRO Astronomy and Space Science and Mount Stromlo Observatory, Canberra ACT 0200 

Australia 

Barney J. Rickett 

Department of Electrical and Computer Engineering, University of California, San Diego, La 

Jolla, CA 92093, USA 

Michael S. Dutka 

The Catholic University of America, 620 Michigan Ave., N.E., Washington DC 20064, U.S.A. 

Jun Yi Koay 

ICRAR, Curtin University, Bentley, WA 6845, Australia 
James E. J. Lovell 

School of Mathematics and Physics, University of Tasmania, TAS 7001, Australia 

Hay ley E. Bignall 

ICRAR, Curtin University, Bentley, WA 6845, Australia 
Lucyna Kedziora-Chudczer 

: Now at NASA/Goddard Space Flight Center/ORAU, Code 661, Astroparticle Physics Lab., Greenbelt, MD 
20771, U.S.A. 

2 Now at Faculty of Engineering and Science, Universiti Tunku Abdul Rahman, Setapak, KL 53300, Malaysia 



- 2 - 


School of Physics and Astrophysics, UNSW, Sydney NSW 2052 Australia 

and 

Jean-Pierre Macquart 

ICRAR, Curtin University, Bentley, WA 6845, Australia 

ABSTRACT 

Intraday variability (IDV) of the radio emission from active galactic nuclei is now 
known to be predominantly due to interstellar scintillation (ISS). The MASIV (The 
Microarcsecond Scintillation Induced Variability) survey of 443 flat spectrum sources 
revealed that the IDV is related to the radio flux density and redshift. A study of the 
physical properties of these sources has been severely handicapped by the absence of 
reliable redshift measurements for many of these objects. This paper presents 79 new 
redshifts and a critical evaluation of 233 redshifts obtained from the literature. We 
classify spectroscopic identifications based on emission line properties, finding that 78% 
of the sources have broad emission lines and are mainly FSRQs. About 16% are weak 
lined objects, chiefly BL Lacs, and the remaining 6% are narrow line objects. The gross 
properties (redshift, spectroscopic class) of the MASIV sample are similar to those of 
other blazar surveys. However, the extreme compactness implied by ISS favors FSRQs 
and BL Lacs in the MASIV sample as these are the most compact object classes. We 
confirm that the level of IDV depends on the 5 GHz flux density for all optical spectral 
types. We find that BL Lac objects tend to be more variable than broad line quasars. 
The level of ISS decreases substantially above a redshift of about two. The decrease is 
found to be generally consistent with ISS expected for beamed emission from a jet that 
is limited to a fixed maximum brightness temperature in the source rest frame. 

Subject headings: galaxies: distances and redshifts — galaxies: intergalactic medium — 
galaxies: ISM — galaxies: BL Lacertae objects : general — galaxies: quasars: general 


1. INTRODUCTION 

The discovery of centimeter wavelength, short term intra-day variability in some compact, 
flat-spectrum extragalactic radio sources (Heeschen 1984; Heeschen et al. 1987) was of immediate 
and profound astrophysical consequence. This phenomenon, which was soon dubbed Intraday 
Variability (IDV; Wagner & Witzel 1995), implied brightness temperatures up to 10 21 K in extreme 
cases (Kedziora-Chudczer et al. 1997) if the variations were intrinsic. Such enormous brightness 
temperatures present a serious challenge to the current paradigm for the physics of extragalactic 
radio sources. 
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However, there is now overwhelming evidence that IDV results primarily from scintillation in 
the turbulent ionized interstellar medium of our Galaxy, commonly referred to as interstellar scin- 
tillation (ISS) 1 . This explanation was proposed by Heeschen & Rickett (1987) and has now been 
confirmed from several further lines of observational evidence. Time delays have been detected 
between the arrival times of the intensity fluctuations from ISS sources at two widely-spaced tele- 
scopes (Jauncey et al. 2000; Dennett-Thorpe & de Bruyn 2002; Bignall et al. 2006). Further, annual 
cycles have been detected in the timescale of ISS sources (e.g. Jauncey &; Macquart 2001; Rickett 
et al. 2001; Bignall et al. 2003; Dennett-Thorpe &: de Bruyn 2003), a periodic modulation that 
most plausibly results from the change in relative velocity of the Earth and the scattering medium 
through the course of a year. Subsequently, a strong correlation has been found between the level 
of IDV and Galactic emission measure estimated from Ha measurements (Lovell et al. 2008). This 
shows ISS to be the predominant cause of radio IDV. Though some questions remain (Jauncey 
et al. 2001; Krichbaum et al. 2002), interstellar scintillation is the only reasonable explanation of 
these observations (Jauncey et al. 2001). 

Though an extrinsic (as opposed to intrinsic) origin for IDV points to less extreme physical con- 
ditions for the sources that exhibit this phenomenon, scintillators are still among the most extreme 
and active radio AGN known. For a source to scintillate, its angular size must be comparable to 
that of the first Fresnel zone (Narayan 1992) which implies microarcsecond angular sizes for screen 
distances of tens to hundreds of parsecs. Such a high resolution cannot be achieved by any other 
existing technique, including space VLBI. Thus, interstellar scintillation has the potential to probe 
within a few light months of the central black hole (Bignall et al. 2003). Further, inferred brightness 
temperatures of a few scintillators are well in excess of 10 14 K (Macquart et al. 2000) which implies 
Doppler factors of several hundred or more (Readhead 1994) and is significantly higher than seen 
in VLBI surveys (e.g. Lister et al. 2009; Ojha et al. 2010). While Doppler factors very similar 
to those found in VLBI are implied from the ISS analysis of the Green Bank radio monitoring 
program of 146 compact radio sources (Rickett et al. 2006), an investigation of the properties of 
AGN that exhibit ISS is of considerable astrophysical interest, both because their extreme physical 
properties make them ideal probes of existing models of AGN physics and they probe phenomena 
at unprecedented resolution. 

The Microarcsecond Scintillation-Induced Variability (MASIV) survey was undertaken by 
Lovell et al. (2003, 2008) in order to construct a large sample of scintillating extragalactic sources. 
A key discovery from the MASIV survey is that both the number of sources showing ISS and their 
level of ISS variations appear to decrease with increasing redshift (Lovell et al. 2008). This result 
was obtained for the 275 (of 443) sources with measured redshifts. While 206 of these redshifts 
were obtained from the published literature, 69 of them were based on preliminary analysis of our 
unpublished optical observations. A primary goal of the present paper is to publish the full analysis 
of those optical observations and an additional 22 sources subsequently observed. 


Tn this paper we call the phenomenon of rapid variable flux density IDV and its cause, scintillation, ISS. 
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The apparent decrease in ISS with redshift is a new cosmologically important result. Its 
interpretation must include the cosmological prediction that the angular size of an object should 
initially decrease with redshift but start to increase above redshifts of about one. Observational 
evidence for this phenomenon has long been sought, but the best evidence to date is still only 
marginal (Gurvits et al. 1999). The MASIV sources were selected for their flat radio spectra and 
many are quasars. Thus any detailed interpretation must also address possible evolution in the 
beamed AGN jets presumably responsible for their very compact diameters. Alternatively, there 
could be scattering effects in the inter-galactic medium that could cause radiation from the more 
distant AGN to be scatter broadened so that their angular sizes are too great for them to scintillate 
in the interstellar medium of our Galaxy. 

A follow-up to the MASIV Survey has been conducted in which the ISS of a subsample of 140 
MASIV sources (70 with z < 2 and another 70 with z > 2) were observed simultaneously at 5 GHz 
and 8.4 GHz over a duration of 11 days using the VLA (Koay et al. 2011). These observations 
provided a means of determining the origin of the redshift dependence of ISS, by examining how 
the effect scales with frequency. This exploited the fact that intrinsic source size effects and scatter 
broadening both have different frequency dependences. The analyses and results were presented in 
Koay et al. (2012). 

Our investigation of scintillators has been greatly handicapped by the lack of source redshifts. 
The MASIV survey found more scintillators among the weaker radio source sample, which also 
showed more extreme scintillation behaviour than their higher flux density counterparts. Since 
redshifts have been measured predominantly for the stronger radio sources, we particularly need 
to measure redshifts for the weaker sources. Redshifts are essential to determine physical proper- 
ties including linear sizes, accurate brightness temperatures and luminosity distributions. Optical 
luminosities are also an important probe of the physics of Compton scattering. For high Lorentz 
factor jets, energy losses should be dominated by Compton scattering of the diffuse radiation field 
of the cosmic microwave background, starlight and reprocessed emission from the nucleus (Begel- 
man et al. 1994). This energy will escape mainly at X-ray wavelengths but some is expected in the 
optical band. Thus redshifts are needed to calculate the luminosities in order to understand the 
total energy budgets of these extreme sources. 

In this paper we present new redshifts obtained with the 2.56 m Nordic Optical Telescope 
(NOT) on La Palma 2 and the 5 m Hale Telescope at Mount Palomar. We also present redshifts 
obtained from the literature whose reliability could be ascertained and optical magnitudes collated 
from the surveys and literature. Further, we look for relationships between the redshift, spectral 
classification, luminosity and the observed level of ISS variation and compare the results with a 
simple model for ISS. 


J Based on observations made with the Nordic Optical Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 



- 5 - 


Throughout this paper we adopt a cosmology with fi m = 0.27, Ha = 0.73 and Ho = 70 km s 1 
Mpc -1 . We use the convention S u oc v a for flux density at frequency v. 


2. THE MASIV SAMPLE 

The MASIV sample has its roots in two radio catalogues: CLASS (Cosmic Lens All Sky Survey; 
Myers et al. 2003) and JVAS (Jodrell Bank VLA Astrometric Survey; Patnaik et al. 1992; Browne 
et al. 1998; Wilkinson et al. 1998). Both of these were targeted surveys of flat-spectrum radio 
sources using the VLA at 8.4 GHz in the ‘A’ configuration. JVAS is effectively a bright sub-sample 
of the complete JVAS/CLASS sample, which consists of more than 11,000 sources with flux density 
down to 30 mJy (Myers et al. 2003). We note that the CLASS flux density limit is well below the 
100 mJy flux density limit chosen for the MASIV sample (see below). 

For MASIV observations, we wanted to target compact flat-spectrum sources that were un- 
resolved at 8.4 GHz with the VLA and located north of the equator. As a first step, all CLASS 
sources with modeled source sizes < 50 mas and all JVAS sources with > 95% of their flux density 
in an unresolved component were selected. These selection criteria were chosen due to the way the 
data were presented in each catalogue: for JVAS, the peak and total VLA flux densities are listed; 
for CLASS the modeled component size is presented. Both criteria select those sources which are 
essentially unresolved point sources for short (0.5 - 1.5 minute) snapshot observations with the VLA 
in ‘A’ configuration (this is the highest resolution configuration of the VLA). Sources where the 
catalogues indicated the presence of any confusing sources - indeed, any emission outside of the un- 
resolved target source - in the VLA field, were dropped from the sample, as selecting each MASIV 
target field to consist of a single, isolated point source greatly simplifies the analysis of variability. 
To obtain a flat-spectrum sample with spectral index a > —0.3 (S„ oc i/ a ), these sources were 
cross-correlated with the NVSS catalogue (NRAO VLA Sky Survey; Condon et al. 1998), although 
we note that the observed flux densities in these catalogues were non-simultaneous. 

The sample was next divided into strong and weak subsamples of about 300 sources each with 
the strong sample consisting of sources with S^aghz > 0.6 Jy and the weak sample consisting of 
sources with 0.103 < S^aghz < 0.13 Jy. Finally, the sources were selected to have uniform RA - 
5 distribution in order to have the best possible coverage for the VLA observations (for details see 
Lovell et al. 2003). 

MASIV’s 5 GHz observations were carried out on this ‘core’ sample of 578 sources. Of these, 
102 sources were removed in the analysis stage due to the presence of structure or confusion while 
one was removed due to an error in the initial sample selection process. This left a sample of 475 
point sources common to all four epochs of MASIV observations. However, as described in Lovell 
et al. (2008), 32 of these sources were used as secondary calibrators in two or more epochs, which 
effectively removes their low-level (instrumental or real) variations in those epochs, relative to the 
other sources. These 32 sources were excluded from the final MASIV sample, leaving 443 sources 



- 6 - 


for full analysis. 

As the average flux densities of many sources from the originally defined radio “strong” and 
“weak” samples vary with respect to the original catalogued flux densities, we used the first epoch 
VLA data (Lovell et al. 2003) to set the dividing flux density to 0.3 Jy at 5 GHz. In the rest of the 
paper we refer to the radio strong sample of 229 sources (including the 32 secondary calibration 
sources for the optical analysis only, Section 4) and the radio weak sample of 246 sources. 

Lovell et al. (2008) report which of these sources exhibited Intra Day Variability (IDV) and 
in how many epochs of observation were classified as variable. A structure function (SF) analysis 
of the flux density variations was also done, which provides a quantitative measure of the strength 
of the variation. The flux density was normalized by its mean and the average SF was computed 
from the 4 epochs and tabulated for each source as D(2days), the SF at a time lag of two days. 
This quantifies the variations and was shown to be correlated with the interstellar electron column 
density, estimated from Ha emission on that line of sight. The correlation confirms that the 
variations are indeed predominantly due to ISS. There are thus two ways of quantifying the ISS 
of a source, either as the number of epochs in which it was classified as variable or the 4-epoch 
mean D(2days), where we note that no significance should be attached to values smaller than about 
4 x 10 -4 (for details, see Lovell et al. 2008). 


3. REVIEW OF THE LITERATURE 

The main goal of this work is to present the spectroscopic identifications and redshifts of the 
flat-spectrum, extragalactic radio sources that make up the MASIV sample. 


3.1. Redshifts from the Literature 

We searched for redshifts and spectroscopic identifications of MASIV objects in the literature 
starting from data repositories like the NASA Extragalactic Database (NED 3 ). Spectroscopic iden- 
tifications and redshifts were found for the majority (171 of 229) of the radio strong sample sources. 
However, this is not the case with the weak sample (62 of 246). Most literature redshifts and spec- 
troscopic identifications are from the 12th Veron-Cetty & Veron (2006) compilation of known AGN, 
hereafter VCV12. Identifications for 27 objects are from CGRaBS (Candidate Gamma- Ray Blazar 
Survey; Healey et al. 2008) and the rest from sources found using NED. 

For those objects which had reported redshifts, we tried to locate their spectra and/or infor- 
mation about their emission lines (line identifications, equivalent widths) from which the redshift, 
estimates were obtained. We classified the sources into categories A, B or C, based on the reliability 


3 http:/ /nedwww. ipac.caltech.edu/ 
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of their redshift estimates. Sources with redshifts obtained from multiple reported emission lines or 
with spectra clearly showing multiple emission lines were placed in category A. If only a single line 
was reported with wide wavelength coverage, or the spectrum provided was of marginal quality, 
the source was placed in category B (see comments in Table 2). Redshift estimates of sources for 
which the spectra were not available, but where the process by which the redshifts were obtained 
were clearly described, were also placed in category B, as were three sources where imaging studies 
have detected the host galaxy but show featureless spectra. For these three objects, we adopted the 
imaging redshift which is estimated assuming the host galaxy is a “standard candle” (for details see 
Sbarufatti et al. 2005). Finally, objects in which the redshifts were indicated to be photometric, or 
where no further information could be found on how they were obtained, were placed in category 
C. Objects with (nearly) featureless spectra, with conflicting reported redshifts, or with redshifts 
estimated based on an intervening absorption system, were all placed in category C. Comments on 
archival redshifts that either were rejected or judged less reliable are in Appendix A. 

3.2. Optical identification 

In order to define the sample for optical spectroscopy and select appropriate observing re- 
sources, accurate optical identification was a prerequisite. 

For optical identification we used the CDS client database 4 to access three large area surveys, 
the Sloan Digital Sky Survey (SDSS), the Guide Star Catalog 2.3 (GSC) and the USNO-B1 Catalog. 
The SDSS DR5 (Adelman-McCarthy et al. 2007) is the primary catalogue for optical identification. 
It is complete down to r=22.5 (A e ff ~6230A) with astrometric accuracy of about Of! I (Pier et al. 
2003) and covers essentially the Galactic gap (RA ~12 ft 6 ~30°) and some smaller patches altogether 
about 8000 square degrees. As a secondary catalogue we used the GSC 2.3 (Lasker et al. 2008) 
which is compiled from scans of the POSS-II plates and has F magnitudes (X e ff 6500A) down to 
~20.5. This has poorer photometric (0.13-0.22 mag) and astrometric accuracy (0 / . / 3-0 / . / 4) than the 
SDSS. For 16 sources which were not found from the GSC2.3 we used USNO-B1 R2-magnitudes 
(Monet et al. 2003), which has poorer photometric accuracy than the GSC2.3 (Sesar et al. 2006). 
Of these sources ten are faint (R2 >19.0), just visible from the POSS-I O-plates, and four are 
bright extended targets (R2 <11.3). The ten faint objects might be missing from the GSC2.3 due 
to variability, while the four bright, extended objects may be missing due to the design of the 
GSC2.3. Finally, for eight faint sources we found optical magnitudes from the literature. 

The initial search radius around the radio position was 5". From SDSS DR5 we found 165 
matches to a MASIV radio source within one arcsecond, with median r-magnitude of 18.7T1.8 
(standard deviation of the distribution) and the faintest optical identifications having r=22.9. The 
mean separation increases from 0707 for r < 20.5 mag sources to 0715 for the 37 sources with 20.5 

4 http:/ /cds. u-strasbg.fr/ 



< r < 22.9. The GSC2.3 has 356 matches within F/5 with median magnitude of 18.44=1.5. The 
mean separation of F < 19.5 sources is 0725 and for the 38 sources with 19.5 < F < 20.7 it is 
0 / . / 39. Fewer than 10% of the GSC2.3 matched sources have separations greater than 0.5 arcsecond 
and only ~2% have separations greater than one second of arc. In Figure 1 we show a plot of 
the measured separation versus apparent magnitude. In both the SDSS DR5 and GSC2.3 matched 
sources, the accuracy of the astrometry is as expected and decreases close to the limiting magnitude 
of the survey. 

Multiple matches were found in 14 cases (~4% of the matches). Inspection of the data and 
visual inspection of the images (available from the CDS) revealed that close companions, defined 
as those with a separation less than one arcsecond (8 sources), were from separate epochs. This 
suggests that there is only a single source, but the astrometric solution is slightly different from 
epoch to epoch. For the remaining six sources the nearer object was selected as the optical coun- 
terpart and a second target was found with the angular distance 2. 5-4. 9 arcseconds from the radio 
source. Of the 356 GSC2.3 detected objects, 133 also have SDSS DR5 magnitudes, hence we are 
left with 223 sources with only F-magnitudes. Down to a Galactic latitude of 5 degrees, all SDSS 
and GSC2.3 matches had only a single source within the search radius. Table 1 summarizes the 
SDSS and GSC identifications. 

In addition we defined a subsample, (08 ft <RA< 16^ and <5 < 64°), which covers roughly 
the SDSS DR5 sky area (hereafter ‘SDSS sample’). Of the 79 radio strong sources, all but one 
secondary calibrator source have optical identifications. Optical identifications are available for 82 
out of 95 sources in the radio weak sample. 


3.2.1. Combining GSC2.3 and SDSS DR5 Data 

We transformed the GSC2.3 and SDSS DR5 data to the common magnitude system using 

Rsdss = r - 0.272(r — i) - 0.159 (1) 

by Chonis & Gaskell (2008) for the SDSS data. This is derived using SDSS DR5 and Landolt 
(1992) data for r > 14 magnitude stars. For the r < 14 stars they found a systematic magnitude 
difference, however we have only two such bright targets. 

We obtained the GSC2.3 F- and N-magnitudes for the standard stars, which were observed 
more than four times by Landolt (2009) where we found (almost) a one to one correspondence 

Rgsc = F. (2) 

The F — N color term and zero point corrections are small and are omitted (see Appendix B for 
details). The USNO-B1 R2-magnitudes were converted in a similar way as the F magnitudes. For 
249 sources with 13 < F < 19 the average magnitude difference F — R2 is ~-0.07. Taking into 



- 9 - 


account the small number of R2-magnitudes and the magnitude offset, the possible bias to the 
global properties is negligible. 

We checked the common magnitude system by selecting all targets with SDSS photometry and 
GSC2.3 F magnitude between 16.5 and 19.5, altogether 97 sources. The median r - F magnitude 
difference is 0.3T0.4 and after the Rsdss ~ Rose transformations the median difference is 0.14=0.4 
(Figure 2) suggesting that the systematic difference between the transformed Rsdss an d Rgsc 
magnitudes is small in comparison to the uncontrollables, such as variability of the sources and 
photographic plate sensitivity changes. 

The magnitudes are integrated magnitudes and there has been no attempt to separate the host 
galaxy contribution from the nucleus or correct the emission line fluxes. For optically extended 
sources (all low redshift, z < 0.2 sources as well as Sy2/NLRG-type objects at higher redshift) the 
absolute magnitude includes both the core and the host galaxy. For example, J0057+3021 (NGC 
315) has USNO R magnitude ~11.1 however the core R-magnitude is ~19.8 (assuming V-R=0.3; 
Verdoes Kleijn et al. 2002). Also, the GSC2.3 photometry pipeline was tuned for point sources, 
thus there is a systematic overestimate for bright non-stellar magnitudes (Lasker et al. 2008). 

All the magnitudes were de-reddened using the E(B — V ) values and the relative extinction 
values from Schlegel et al. (1998). In Table 2 and for the optical luminosities, we refer to Rsdss 
and Rqsc as R-magnitude. 


3.2.2. Radio and Optical Properties 

The optical identification rate, combining the results from the SDSS DR5 and GSC2.3 cata- 
logues, is ~ 86% for the radio strong and ~ 77% for the radio weak sample. However, since the 
two catalogues have different sky coverage and hence different Galactic latitude coverage, this could 
bias the identification rate. Therefore, using only the (all sky) GSC2.3 catalogue, the strong sample 
has ~ 80% (187) of F-magnitudes brighter than 20 compared to ~ 60% (146) for the weak sample. 

Including empty fields 5 , the median optical flux density of the strong radio sample is ~0.15 
mJy (18.2 mag) decreasing by a factor of 3.2 to ~0.047 nrJy (19.5 mag) for the weak sample. In 
comparison, the median radio flux density decreases by a factor of 7.4, from ~0.89 Jy to ~0.12 
Jy in the two samples. If we were to use this radio flux ratio to calculate the weak radio sample 
median optical magnitude from the strong radio sample mean optical magnitude, it would yield a 
R-magnitude of 20.1 for the weak sample. The /^-correction difference between the weak and strong 
samples is unknown, but most likely it is ~0.2 magnitudes or less on average (e.g. Wisotzki 2000), 
assuming the average redshifts of 1.2 and 1.8 for the strong and weak samples respectively (see 
below). 


'In the optical there is a well defined upper detection limit, i.e. the median is a non-biased estimator. 
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Figure 3 shows the optical versus radio flux density with diagonal lines indicating constant 
radio-optical spectral index (defined by S v oc ^“ ro ), which has only been corrected for Galactic 
extinction and lacks any ^-correction because 40% of the radio weak sample objects have unknown 
redshift. The top line corresponds to a ro =— 0.2, the generally accepted radio quietness limit (Stocke 
et al. 1985). Although there appears to be no detailed correlation between optical and radio flux 
density, the weak sample is on average one optical magnitude fainter than the strong sample. The 
radio weak sample also appears to have an even larger spread than does the radio strong sample. 
The weak sample has 17 sources with a ro > —0.45 and the strong one only five. Close to ov o =-1.0 
both samples have many optically empty fields thus we cannot draw any conclusions about the 
underlying distribution. It is also important to note that the non-simultaneous nature of these 
radio and optical observations will further increase the scatter. The difference between the radio 
to optical flux ratios of the radio strong and weak samples is statistically significant. Running the 
ASURV Kaplan-Meier estimator (Avni et al. 1980), the probability that the two samples are drawn 
from the same parent population is small (P <0.01%). This holds for the ‘SDSS sample’ and also 
for the whole sample using a shallower GSC2.3 limiting magnitude. 

Our results indicate that the radio faint objects are optically “brighter” than expected using 
the radio strong objects as reference. The explanation for the difference is not clear. One possibility 
is that the radio weak sources have larger host galaxy contributions in optical than radio strong 
objects. Another possibility is that the average spectral energy distributions of the strong and 
weak sources are different, perhaps due to differing contributions of the accretion disk (“big blue 
bump”) to the optical flux density or different amounts of dust. 


4. SPECTROSCOPIC FOLLOW-UP 

The new redshifts reported in this paper were obtained using the 2.56 m Nordic Optical Tele- 
scope (NOT 6 ) which is located at Roque de los Muchachos, La Palma, Canary Islands, Spain and 
the 5 m Hale Telescope at Mount Palomar, California, USA. Most of the NOT data were obtained 
during two observing runs (July 3rd to 7th, 2005 and July 22nd to 28th, 2006) supplemented by 
a few additional nights between January 2004 and July 2007. The Mount Palomar data were ob- 
tained during an observing run from August 9th to 10th, 2007. Most fields were pre-imaged at 
NOT with deeper and better image quality than the SDSS frames to verify the optical identification 
of the radio source. Altogether 90 objects were observed at NOT and 28 objects were observed at 
Palomar with some targets observed at both telescopes. 

At NOT, low resolution spectra were obtained using the Andalucia Faint Object Spectrograph 
and Camera (ALFOSC 7 ) with Grisrn 4 (300 rulings/mm, 3200-9100 A, dispersion of 3A/pixel) with 


6 http:/ /www. not. iac.es 

7 The data presented here have been taken using ALFOSC, which is owned by the Instituto de Astrofisica de 
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second order blocking filter giving wavelength coverage about 3800 -9100 A. However, in practice 
the NOT red limit of the spectrum is reduced to about 8000 A due to strong fringing. We used 
Grisrn 6 (600 rulings/mm, 3200-5550 A) or Grism 3 (400 rulings/mm, 3200-6700 A) for objects with 
a single detected line at about 5000 A in order to confirm the redshifts from either Mgll, CIII] or 
the Ly-alpha line. The NOT run in 2005 was photometric, based on Carlsberg Meridian Telescope 
(CMT) data 8 , with seeing of 0 / /6-l / /4 and the 2006 run had some dusty and/or non-photometric 
nights with seeing. At Hale, the Double Spectrograph (DBSP) was used with a setup of 

600 lines on the blue arm (3400-5400 A, ~lA per pixel) and 158 lines for the red (5400-9500 A, 
^4. 9A per pixel) with either l"A) or T/5 slit giving continuous wavelength coverage from about 
3400 to 9500 A. 

Most of the targets were observed with the slit in parallactic angle in order to reduce the light 
losses in blue. We note that, each spectroscopic target has a characteristic AGN spectrum and there 
is no evidence of any spurious counterparts. In two cases the nearby companion was put on the slit 
and in both cases the object closer to the radio source had a characteristic AGN spectrum while 
the companion had a typical stellar spectrum. For the NOT data, an internal flat field (Halogen) 
image was obtained before each science exposure, in order to improve the fringe correction at the red 
end. In addition, several HeNe images were taken during the night for the wavelength calibration. 
The DBSP calibration frames, dome flats, HeNeAr (red arm), FeAr (blue arm) images, were taken 
morning and evening. The data were reduced using standard IRAF 9 procedures including bias 
subtraction, flat field correction using internal flat field lamp and wavelength calibration. Before 
extracting the spectrum, cosmic rays were removed from individual spectrograms. The majority of 
the spectroscopy targets have R-magnitude between 17 and 21, with a median of 18.8. 

The redshifts were determined from the narrow emission lines (e.g. [O III]AA 4959, 5007, [O II] A 
3727) whenever possible, or from strong broad emission lines otherwise, typically H (3, Mgll , C III], 
C IV or Lya. The new redshifts are based on two or more line identifications. However, in four cases 
with wide spectral coverage and high S/N, we assumed the line to be Mgll , as if the line were to be 
C III], we would have expected to detect C IV] or Mgll . In Table 2 we list, respectively, source name, 
RA, declination, first epoch 5 GHz flux density, epochs of variability (at 5 GHz), D(2days), Galactic 
reddening, de-reddened R-magnitude, R-magnitude reference, redshift, spectroscopic identification 
with the original references and a comment. Selected objects are discussed in Appendix A and C. 
Objects with only a single detected emission line are listed in Table 3. 


Andalucia (IAA) and operated at the Nordic Optical Telescope under agreement between IAA and the NBIfAFG of 
the Astronomical Observatory of Copenhagen. 

8 http:/ /www. ast.cam.ac.uk/~dwe/SRF/camc_extinction.html 

9 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation (NSF). 
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4.1. Classification Scheme 

On the basis of the optical spectrum the sources have been divided into three main groups: 
Type 0 objects with weak emission lines (mainly BL Lacs), Type 1 objects with strong and broad 
emission lines (mainly FSRQ) and Type 2 objects with strong and narrow emission lines (mainly 
Seyfert 2’s). In the context of AGN unification schemes, Type 0 and Type 1 objects have their 
Doppler boosted relativistic jet aligned close to the line of sight with Type 2 objects seen at a 
larger angle in such a way that the broad line region is obscured by a dusty torus (see review by 
Urry & Padovani 1995). The spectroscopic classification scheme is adapted from previous works 
e.g. Laurent-Muehleisen et al. (1998), Caccianiga et al. (2002b) and Caccianiga et al. (2008). Many 
observed properties of the AGN, such as emission line equivalent widths ( EW ) and luminosities 
have continuous distributions. Hence the division between, for instance, Type 0 and Type 1 based 
on equivalent widths EW or between Seyfert 1 and FSRQ based on optical luminosity depend on 
the dividing line chosen. The VCV12 classifies all objects fainter than Mr = —23 (using VCV12 
cosmology and magnitudes) as (low luminosity) AGN, which is converted to our cosmology and 
R-band resulting in a QSO vs radiogalaxy (RG) limit Mr = —22.8. As noted by the VCV12, the 
Mr =-23 is an arbitrary limit without any physical reason. Due to different source classification 
schemes, and/or to variability, a QSO might have been classified as an AGN in one epoch but as 
an FSRQ at another epoch, and vice versa. For example, a strong radio source, J 1506+3730 is 
classified as an AGN by Healey et al. (2008), NLRG by Sowards-Emmerd et al. (2004) and Sy2 by 
VCV12. Also, in some cases the classification depends on S/N, e.g. a weak AGN-core with dilution 
by the host galaxy where broad emission lines might be undetected; for J 0057+3021, Ho et al. 
(1995) reported a broad H a line, but a later study could not confirm it (Gongalves & Serote Roos 
2004). 

We describe our classification scheme below. Note that all the line equivalent widths (EW) 
are rest frame values: 


0a) BL Lacs: objects with no lines or line rest frame EW < 5A and C < 25%, where contrast C 
is defined as 

C = (/+-/_)//+, (3) 

where /_ is the mean flux in the 3750-3950 A range and /+ in the 4050-4250A range where 
flux is measured in frequency. In the case of moderate S/N we give only tentative spectroscopic 
identification. 

0b) BL Lac candidates: 25 < C < 40% and strongest line EW < 40 A. 

0c) Passive elliptical galaxies (PEG): C > 40% and EW < 40A. 

la) Flat spectrum radio quasar (FSRQ): at least one permitted line excluding Ha , line FWHM > 
1000 krn/s and Mr < -22.8. 

lb) Seyfert type 1 galaxies: Mr > -22.8, the subclass has been defined based on the H/3/[ O III] A5007 
line ratio (Winkler 1992). 
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lc) Narrowline Seyfertl/QSO (NLSyl/NLQSO): At low redshift: Hf }/[ O III]A5007 < 3 and Hf3 < 
2000 krn/s. For the z > 1.5 targets: rest frame line width 1000-2000 km/s and the relative strength 
of the He II emission line compared to C IV following Heckman et al. (1995); Caccianiga et al. 
(2008). Type 1 objects have weak He II emission line. 

2a) Narrow line Radio galaxy (NLRG): line FWHM < 1000 km/s and Mr < -22.8. 

2b) Seyfert type 2 galaxies: Mr > -22.8, the subclass has been defined based on the H/3 /[ O III] A5007 
line ratio (Winkler 1992). 

2c) Galaxy: C > 40%, absorption lines e.g. Mg I and no emission lines 

Low redshift sources are the most problematic to classify as host galaxy dilution can hide 
broad lines and can reduce the line EW. Also the non-thermal continuum can be difficult to detect. 
Fortunately only four of our follow-up targets are at low redshift (z < 0.2), suggesting that the bias 
from the uncertain spectroscopic identification is small. 


4.2. Summary of the spectroscopic follow-up 

We report 79 new redshifts with spectroscopic identifications, 7 new BL Lacs with featureless 
spectra, and twelve objects with a single emission line of which six have a tentative redshift (Table 3 
and Appendix C). We confirm the BLLac status of nine objects, and confirm the archival redshifts 
as well as spectroscopic identifications for nine other objects. We also re-classify two BL Lac objects 
using our own data (Appendix A). Spectroscopic identifications from the literature were reviewed 
using the same scheme as for our own data, whenever possible. In some cases this led us to adopt 
a different classification than NED (see appendix A). Example spectra of our own observations are 
shown in Figure 4. 


4.3. RESULTS AND SAMPLE PROPERTIES 

Here we present 347 MASIV sources which have their spectra classified, either from the lit- 
erature or based on our own data (Table 2). The radio strong sample is 91% (208/229) complete 
and the radio weak sample is 57% complete (139/246). Of the twenty-one radio strong sources 
without spectroscopic identification, seven have magnitudes from the surveys or literature, with 
R-magnitude between 17.3 and 23.99. 

In Table 4 we summarize the number of sources in each optical spectral classification, separating 
the radio strong and weak subsamples. It shows that Type 1 sources are the most common, 
comprising 78% of both the radio strong and weak samples. About 18% of the radio strong sources 
and 16% of the radio weak sources have Type 0 spectra. A small minority of the sources have Type 
2 spectra, 4% and 6% of the radio strong and radio weak samples, respectively. It also shows that 
almost all Type 1 sources are optically luminous FSRQs. In Table 5 we compare the identifications 
based on the present sample and the ‘SDSS sample’. The results suggest that the two samples have 
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similar distributions of the sources for both radio strong and weak subsamples. Also the redshift 
distributions are identical. This suggests that the observed sample is representative in comparison 
to the full 475 object sample. 

Though broadly similar, there is some difference in the spectroscopic identifications of the radio 
strong and weak samples. The radio weak sample appears to have slightly more low luminosity 
AGN (e.g. PEG, LINER, Hll-region, Galaxy) and fewer BL Lac objects than the radio strong 
sample. An increasing number of low luminosity objects are expected to be detected as the flux 
limit of the sample decreases. At low redshift (z < 0.2) most (7 out of 8 ) of the radio strong sample 
objects are BL Lacs or have broad emission lines, in contrast to the weak sample, where half (6 
out of 12 ) of the z < 0.2 objects have narrow emission lines or “galaxy” type spectra and a further 
five are classified as PEG. Also, the PEG and galaxy type objects are absent in the radio strong 
sample. 

All together, redshift information is available for 319 sources. Figure 5 shows the redshift 
distributions of subsamples divided by the optical spectral classification and the MASIV first epoch 
5 GHz flux density. For Type 1 sources the median redshift increases from 1.3 to 1.4 as the median 
(5 GHz) flux density decreases from 0.86 to 0.12 Jy. The Type 1 radio strong sample has a clear 
maximum in the distribution, but the radio weak sample has a fairly flat distribution between 0.7 
and 2.0. Of the Type 1 sources ~20% of the radio strong and ~29% for the radio weak sample 
have z > 2.0. The redshift trend due to the flux density limit is similar to that seen from the 
previous blazar surveys e.g. from S5 z mea n= 1-18 (limit 1 Jy) to S4 z m ean= 1-3 (limit 0.5 Jy) to 
DXRBS Zmean= 1-56 (limit 0.05 Jy) (Landt et al. 2001). Most of the Type 0 or Type 2 sources have 
measured redshifts of less than one and almost 50% of the Type 0 sources are without a redshift. 

We note that there might be a small number of mis-classified objects. As in every survey, the 
signal to noise ratio of the spectroscopic data is not always optimized for spectroscopic classifica- 
tion but for measuring the redshift. However, in most cases there is no ambiguity in our source 
classification. As our main results and conclusions are based on Type 1 sources (see below), the 
effect of some mis-classified objects is not significant. 


4.4. MASIV and other Blazar Surveys 

In Table 6 we compare spectroscopic identifications, limiting magnitudes and the radio flux 
density limit of MASIV with some other blazar surveys. The DXRBS sample (Padovani et al. 
2007), is selected by cross-correlating X-ray sources with radio sources (S^ghz <1 50mJy). The 
CGRaBS (Healey et al. 2008) has flat spectrum radio sources selected from the CRATES (Healey 
et al. 2007) survey with S^&ghz > 65 mJy, where the aim is to provide a catalogue of likely 7 -ray 
loud AGN. This selection is based on a figure of merit (FoM) which measures the likelihood that an 
individual radio /X-ray source is associated with a known 7 -ray source. In terms of selection criteria 
the CLASS blazar survey (CBS) (Marcha et al. 2001; Caccianiga et al. 2002b) is most similar to 
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MASIV. The CBS has flat radio spectrum sources with the weakest sources S^ghz ^ 30mJy and 
red magnitude < 17.5. It is interesting to ask how the source identification differs between the 
MASIV selection criteria in comparison to the previous blazar surveys. 

The DXRBS classification into the categories of “Radio Loud Quasar, BL Lac and Radio 
Galaxies” is almost identical to our Type 1, Type 0 Type 2. The possible difference is that they 
might classify “PEG” as “Radio Galaxy” and their survey is nearly complete (~95%, Padovani 
et al. 2007). The CGRaBS classifies objects as FSRQ (Type 1), AGN (Type 1) BLL(Type 0), 
NLRG(Type 2) and GALAXY (Type 2) where our corresponding classifications are in parenthesis. 
CGRaBS is about 79% complete with respect to the entire survey and 85% for objects with known 
R-mag< 23. The comparison between the CBS (~91% complete; Caccianiga & Marcha 2004) and 
MASIV is straightforward, as the MASIV classification is adopted from the CBS. 

The distribution of the spectroscopic classes is almost identical between MASIV and the 
DXRBS and the CGRaBS with slightly more broad line quasars in MASIV than other surveys. 
The main difference between the CBS and MASIV is that the CBS has more low redshift Type 
2 objects than MASIV. This difference remains when using similar radio (‘CBS rb’) and optical 
flux density (‘MASIV ob’) limits for these samples. This indicates that the selection of compact 
radio structure in defining the MASIV sample, filters out nearby Type 2 objects and favors FSRQs. 
Probably as a result of its bright optical flux density limit, the CBS has at least 40% of low redshift 
(z <0.15) sources. In contrast, MASIV has 4% for the total sample, 1% and 6% for the radio strong 
and weak samples, respectively. In all these cases it is assumed that the BLLacs without redshift 
have z >0.15. Another difference between MASIV and the CBS is that sources that are optically 
less luminous but powerful in the radio (Figure 6 top panel, where Mr > —23 and P > 10 26 W/Hz) 
are missing from the optically bright CBS sample. 

The source and the redshift distributions of the MASIV, DXRBS and CGRaBS samples are 
similar. There might be slightly fewer z > 3 sources in CGRaBS (~ 3%) than in MASIV (~5%). 
The difference could be, for example, due to incompleteness of the MASIV sample or pre-selection 
of only compact radio sources for the MASIV sample. MASIV has fewer low radio power FSRQ 
sources (P r < 10 25 " 5 W/Hz) than DXRBS (22%; Landt et al. 2001), but about the same amount 
(~2-3%) as the combined 1-Jy and S4 samples. Some of the sources classified as “galaxy” might 
have weak broad lines which are not detected due to dilution by the host galaxy continuum and/or 
insufficiently high signal-to-noise ratio of the data. 

Apart from decreasing the number of low radio power sources, using compact radio core selec- 
tion increases the number of objects in the MASIV sample that have featureless optical spectrum 
i.e. are classified as BL Lacs. Specifically, it is interesting to note that IDV4 (i.e. IDV detected in 
all four MASIV epochs) is one of the most effective preselections to find BL Lac objects. The use 
of persistent IDV as a selection criterion strongly favours BL Lac objects with 43% of the MASIV 
sources showing IDV in all four epochs identified as BL Lac objects. In comparison, a radio-X-ray 
selection in the XB-REX sample (Caccianiga et al. 2002a) and DXRBS (Padovani et al. 2007) have 
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a BLLac fraction of 25% and 18% respectively. However, 7 -ray loudness may be an even stronger 
selector of BL Lacs with about 50% of Fermi/LAT identifications being BLLacs (Abdo et al. 2009). 

Lister et al. (2001) studied the IDV properties of the Pearson-Readhead compact extragalactic 
radio sources and found that IDV sources have smaller emission line widths and lower 5 GHz 
luminosites than non-IDV sources. Their results about line widths are consistent with our findings 
that 25 to 40% (for IDV1 through IDV4 subsamples) of the MASIV IDV sources are BL Lac 
objects and that 70% of the BL Lacs are IDV sources. Whether there is a correlation between 
IDV and line width amongst the Type 1 sources will be studied in a future paper. The 5 GHz 
luminosity difference found by Lister et al. could be due to combining FSRQs and BL Lacs into one 
IDV/non-IDV sample, and the different redshift distributions of the IDV and non-IDV samples. 


5. Discussion 

5.1. Dependence of ISS on 5 GHz flux density and optical spectral type 

The correlation between IDV and the 5 GHz flux density was discovered in the MASIV survey 
(Lovell et al. 2003) and later confirmed by Lovell et al. (2008). In the latter work, the variability 
of the flux density was studied using two complementary methods (see Section 2): D(2days) which 
combines all the data in a statistical estimator; and a visual analysis of each light-curve at each 
epoch, which classified each source at each epoch as variable or not variable and so allowed an 
estimation of the duty cycle in the variability. Sources were classified as “ISS” if they were found 
to be variable in 2 or more epochs and as “non-ISS” otherwise. 

To study the flux density/ISS correlation in detail, we divided the sample by the flux density 
and spectroscopic type. Figure 7 shows the D(2days) histograms of the MASIV subsamples. There 
is a clear increase of the median and mean D( 2 days) as the flux density decreases and this is seen in 
all spectral types separately. One can see that many of the distributions in Figure 7 are significantly 
skewed, having a tail extending to high values. When binned on a linear scale the distributions 
are more strongly skewed. A measure of this skewness is that the mean values of D(2days) are 
considerably larger than the median, since the rare high values influence the mean but have little 
effect on the median. 

The Kruskal- Wallis (KW) test indicates a ~0.4% probability that the radio strong and radio 
weak Type 0 samples have equal median D(2days). The probability that Type 1 radio strong and 
weak samples have similar median is negligible ( Pkw < 0.01%). Similarly, the Kolmogorov-Smirnov 
(KS) test indicates that the Type 0 radio strong and weak distributions of D(2days) are drawn 
from the same parent population with a probability of 0.1%; for Type 1 sources this becomes pxs < 
0.01%. For Type 2 objects the KW and KS tests indicate the same median or parent population 
with less than 5% probability, however the full sample contains only 19 objects and their D(2days) 
values are typically at or below the threshold (0.0004) for significant variation. This indicates that 
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the D{ 2days) has a S^ghz dependence and this holds for Type 0 and Type 1 sources with high 
statistical confidence. A similar, but less significant, trend is also seen for Type 2 sources. We note 
that we have similar results when comparing the 5 GHz flux density distributions for sources drawn 
from the ‘SDSS sample’. Thus we confirm the IDV-S^G/fz flux density dependence found earlier 
and show that this correlation is present at least for Type 0 and Type 1 sources. 

In previous studies of IDV, the optical spectroscopic type of the objects has not been consid- 
ered. Figure 7 and Table 7 suggest that the IDV properties depend not only on the 5 GHz flux 
density, but also on the optical spectral type. Although the statistical tests indicate a clear differ- 
ence between the distribution of D(2days) in the Type 0 and Type 1 sources for the full sample, 
there are interdependencies that could cause an apparent correlation. 

To study the apparent correlation between spectral type and ISS, we first ask if the different 
redshift distributions of Type 0 and Type 1 sources could cause the difference in the D(2days). 
From previous studies, BL Lacs and FSRQs have been found to have different observed redshift 
distributions (e.g. Giomrni et al. 2012) with BL Lacs having a lower median redshift than FSRQs. 
As many (25 out of 59) of our Type 0 sources are without known redshifts, it is probable that 
the redshift distributions of the Type 0 and Type 1 sources are different. Hence it is difficult to 
compare Type 0 and Type 1 sources accurately. 

First we compared the D(2days) distributions of Type 0 sources with known redshift (z max ~ 
1.15, N=16) to Type 1 sources with ^ <1.15 (N=57). The median D(2days) is higher (0.00127 
vs 0.00041) for the radio strong Type 0 sources than Type 1. The statistical tests indicate some 
difference between the samples; pxs = 0.4% that the distributions are similar, and pkw = 2.5% 
that the medians are the same. When the Type 1 redshift cutoff is increased, the difference between 
the two samples increases. No comparisons were made using sources only from the SDSS sample, 
as there are only nine Type 0 sources; however the median D(2days) of Type 0 sources is higher 
than that of Type 1 sources. 

Next we compare all Type 0 radio strong sources (N=31, median D(2days)=0. 00070) to radio 
strong Type 1 sources with different redshift cutoffs. The KS test indicates p = 2.3% probability 
that the Type 1 with z <1.15, and Type 0 D(2days) distributions are drawn from the same parent 
population and the KW probability that the medians are the same is pkw = 2.3%. When the Type 
1 redshift cutoff is increased, the KW and KS probabilities that the medians or distributions are 
similar decreases. For example, comparing all Type 0 radio strong sources to radio strong Type 1 
sources with z <1.3 the Pkw = 1-1% and pxs = 0.5%. Considering only radio strong sources from 
the SDSS sub-sample, the median D(2days) of the Type 0 is also higher than that of the Type 1 
sources, but the difference is not statistically significant for these smaller samples. 

For the radio weak samples, including all Type 0 and Type 1 sources, there is a statistically 
significant difference between the two samples (pks = 0.3% pkw = 1.8%). However, including 
only z < 1.5 Type 1 sources the distributions are similar. Finally, excluding the low redshift (z < 
0.16) Type 0 sources, mainly PEG, and including only z < 1.5 Type 1 sources the statistical tests 
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indicate significant difference between the two samples (pks = 0.3% pkw = 2.1%). 

Using the current data, we conclude that the Type 0 sources have higher D(2days) values than 
Type 1 sources. However, the statistical significance is dependent on the sample selection, especially 
the redshift cutoff of Type 1 sources. As many Type 0 sources are without known redshift, it is 
impossible to draw solid conclusions regarding whether there is a real D(2days) difference between 
the Type 0 and Type 1 sources. In general, the redshift distributions of BLL (Type 0) and FSRQ 
(Type 1) are different (see e.g. Giommi et al. 2012), with BL Lacs having a lower median redshift 
than FSRQs. Giommi et al. (2012) proposed that BLL and FSRQ would have similar intrinsic 
redshift distributions, however, in this case the D(2days) distributions of Type 0 and Type 1 do 
not match, suggesting different D(2days) properties between Type 0 and Type 1 sources. 

Table 8 shows the number of sources of each spectral type compared with the number of 
MASIV Survey epochs in which the source varied according to the visual classification presented 
by Lovell et al. (2008). This suggests that Type 0 sources show ISS more often than the Type 1 
sources, which is consistent with the D(2days) results above. It is notable that Type 0 sources have 
the highest levels of ISS in both the weak and strong groups, which implies that they typically have 
more compact radio structure than Type 1 sources. The ISS intermittency is likely to be source 
related, rather than purely due to ISM intermittency. These results are relevant to attempts to 
understand the physical differences between the two groups. 


5.2. ISS versus Redshift 

Lovell et al. (2008) found a redshift dependence in the ISS properties. The original 275 radio- 
selected MASIV sources with redshifts revealed that ISS of compact sources at 5 GHz decreases 
with redshift above ~2. From the optical data presented here we can examine whether this decrease 
depends on spectral classification and re-examine the decrease with redshift using the larger sample 
of 320 objects with redshifts. As the D(2days) correlates with the 5 GHz flux density, if the 
5 GHz flux density were to correlate with redshift, this could cause an artificial D(2days)-redshift, 
correlation. We ran Pearson and Spearman correlation coefficient tests for the full and ‘SDSS 
sample’ radio strong and weak Type 0 and Type 1 sources. Out of eight correlation coefficients 
only one was statistically significant, the full radio strong sample Pearson test (r ~ -0.18; p ~3%). 
All the other tests indicated no correlation between S^ghz and redshift. 

The 2008 MASIV analysis of D{ 2days) versus redshift included all sources regardless of their 
spectral classification, specifically including Type 0 sources with known low redshift and excluding 
the Type 0 sources with unknown, possibly high, redshift. We ask whether the presence of Type 
0 sources could be partly responsible for the drop in ISS with increasing redshift. We found (see 
above) the Type 0 sources to be among the strongest scintillators and also to be concentrated at 
redshifts less than 0.7, however almost 50% (27) of the sources do not have a redshift. Fortunately, 
the redshifts of the Type 0 sources are less than 1.15, so the possibly misclassified Type 1 sources 
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will most likely affect the low-redshift end of the Type 1 sources. Lovell et al. (2008) show that the 
average D(2days) values drop around z ~ 2, hence the misclassified Type 1 sources should have 
negligible bias to the ISS-redshift correlation. 

We now restrict the analysis to the Type 1 objects which constitute the great majority of 
our sources. Figure 8 shows the D(2days) against redshift, with horizontal lines indicating the 
median D(2days) for different redshifts bins. Both radio strong and weak sources have the lowest 
.D(2days) median at the highest redshift bin (z > 2.5). Combining the radio strong and weak 
samples, the KW test indicates 0.5% probability that the 1.5 < z < 2.5 and z > 2.5 samples have 
similar D(2days) median and the KS-probability that the two samples are drawn from the same 
parent population is 3.2%. We also note that the two highest redshift bins (H(2days) = 0.00029 and 
0.00051) have about equal numbers of radio strong and radio weak sources, reducing the probability 
that the H(2days) decrease with redshift is induced by the 5 GHz flux density. Figure 8 suggests 
that the radio weak sample has a stronger redshift dependence than the radio strong sample. The 
KW probability that the weak sample H(2days) medians at z > 2.5 and 1.5 < z < 2.5 are the same 
is 2.0% and the KS for the same parent population is 1.4%. From these results we conclude that 
Type 1 sources show a decrease in D(2days) with increasing redshift, though the trend is not as 
strong statistically as the decrease versus 5 GHz flux density. 


5.3. ISS and Radio Power 

We searched for correlations between radio power and the H(2days), and compared the radio 
powers of ISS and non-ISS sources. We calculated the rest-frame radio power at 5 GHz for isotropic 
emission as 

Pr = 4 ttD 2 l S 5 (1 + z)-( 1+a r\ (4) 

where Dl is the luminosity distance in the adopted ACDM cosmology, S$ is the mean first epoch 
MASIV VLA flux density and the last term is for the ^-correction. We have assumed a r =0 for all 
our sources. 

As the initial sample selection has a gap in the 5 GHz flux density distribution between about 
0.2-0. 5 Jy, the radio power distribution is not continuous at a given redshift. In order to study 
continuous distributions, we selected Type 1 sources with 1.0 < z < 1.8, 59 radio strong and 39 
radio weak sources, and divided the sample into ISS sources (with Z4(2days) > 4 x 10~ 4 ) and non- 
ISS sources (D(2days) < 4 x 10~ 4 ) to compare the radio power distributions. The KS and KW tests 
suggest significant difference in the radio power distributions (pks = 0.9 %,pkw = 0.3%), non-ISS 
sources being more powerful. The Pearson correlation coefficient indicates weak correlation (r=- 
0.28 p= 0.5%) between radio power and D(2days). However, as there are fewer radio weak sources 
than radio strong sources (Table 9, Figure 6), this can cause a correlation between ISS and radio 
power. Thus we divided the sample into radio strong and radio weak and compared the radio power 
distributions of the ISS and non-ISS sources. The KS and KW tests show no difference between 
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the samples, suggesting that the correlation seen by the full sample is weak or non-existent. The 
redshift distributions of the samples are identical in each case and the results are not sensitive to 
the selected ISS/non-ISS limit. Our interpretation of the results is that the ISS-P r correlation is 
weak. 

Finally, we note that using the radio loudness limit of P r = 10 23 ' ' W/Hz, (similar to Laurent- 
Muehleisen et al. 1998) five low redshift (z <0.05) sources are “radio quiet” 10 . 


5.4. ISS and Optical Brightness 

In the optical, we searched for correlations between H(2days) and optical magnitude/optical 
luminosity. We selected subsamples, as above, based on the radio power and carried out the 
Spearman and Pearson correlation tests. Apart from Type 0 radio weak sources we found no 
evidence of correlation between the optical R-magnitude and H(2days). As several Type 0 sources 
have low redshifts, and might have a strong host galaxy contribution to the optical magnitude, we 
excluded the low redshift sources (0.01 < z < 0.15), and the indication of correlation disappeared. 
We calculated the rest-frame optical power as in equation 4 but used the R-band flux density and a 
^-correction — 2.5(l+a 0 )logio(l+z), where we assumed a 0 =— 0.5 for all our sources. The statistical 
tests (KS, KW, Spearman) indicate no difference in the optical luminosities of the ISS and non-ISS 
sources nor any correlation between the D(2days) and optical luminosity. 


5.5. Modeling the ISS dependence on Redshift 

The 275 radio-selected MASIV sources with redshifts known prior to the present paper revealed 
that ISS of compact sources at 5 GHz decreases with redshift above z ~ 2 (Lovell et al. 2008). In 
Figure 3 of Koay et al. (2012) we considered the possible influence of a cutoff in radio power or in 
redshift and found that either a single cutoff beyond redshift z=2 explains the data, or different 
cutoffs in P r would be required for the weak and strong samples. Here we follow the simpler 
hypothesis of a physical model in which H(2days) is governed by S^ghz and redshift regardless of 
P r . The model is motivated by the physics of jet emission as discussed below and also by Koay et al. 
(2012). This model also explains why the radio weak sample has a stronger redshift dependence of 
H(2days) compared to the radio strong sample (see Figure 9 in Koay et al. 2012). 

In Figure 9 we have averaged H(2days) into four redshift bins for the FSRQs alone, combining 
both weak and strong sources. We note that excluding the Type 2 sources increases the level of ISS 
in the lowest redshift bin and excluding the Type 0 sources decreases the level of ISS in the lowest 
redshift bin, when compared with our earlier results (Lovell et al. 2008). We also over plot the 


10 J0248+0434, J1141+5953 and J1604+1744 are classified as galaxies, J1112+3527 as PEG and J0057+3021 as a 
LINER 
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data with curves from a simple theoretical model. The model makes many assumptions and so is 
quite simplistic; however, it has the merit of including the most obvious cosmologically important 
parameters. We assume that the scattering is due to a single layer of isotropic Kolmogorov turbu- 
lence in the interstellar electron density centered at 500 pc from the Earth. We further assume that 
the scintillation for a point source at 5 GHz is at the transition from weak to strong scintillation. 
We model the compact part of the sources by a relativistic jet pointed toward the Earth emitting 
beamed radiation characterized by a brightness temperature in the jet frame of T^o = 10 11 K which 
delivers an observable flux density of 1 Jy. This provides a crude model for synchrotron radiation 
that is self-absorbed causing an upper limit in brightness. The jet is characterized by a Doppler 
factor r which is assumed independent of z. The apparent brightness temperature of the source at 
redshift z is thus given by Tj,oT/(l + z). 

Since for a flux density limited sample the apparent angular diameter varies as the inverse 
square root of the apparent brightness temperature, the (1+z) reduction factor causes the apparent 
diameter of the compact source to increase as \/l + z, which corresponds to the cosmological “an- 
gular size versus redshift” phenomenon. VLBI observers have long sought to detect this predicted 
increase in angular diameter with redshift, but with only partial success. For example, Gurvits 
et al. (1999) show angular diameters of compact radio sources that exhibit a barely significant 
increase beyond a redshift of about 1. We see a highly significant drop in D(2days) with redshift, 
but a model for the ISS phenomenon is required before it can be interpreted cosmologically. An 
ISS model was applied to the data from the Green Bank monitoring observations by Rickett et al. 
(2006). However, here we have improved the model by using the formulae developed by Goodman 
(2006), which explicitly include contributions from diffractive scintillations, that are important near 
the transition from weak to strong scintillation and were not included by Rickett et al. (2006). An 
understanding of D(2days) is obtained by noting its asymptotic dependence oc 0~^ irce for intraday 
ISS. This assumes that the ISS is quenched and that its time scale is not much longer than 2 days. 
However, if the time scales are much longer, then D(2days) drops even more steeply with Source. 
The values chosen for T),o, T and the mean flux density S c in the compact source enter through a 
single quantity T^T / S c . Thus, while the model exhibits the expected dependence on redshift, we 
should also expect a wide range in this unknown quantity. 

A single model curve (T = 20) can be seen to overlap three of the four observed values within 
Tier. However, in the highest redshift bin, D(2days) lies more than 2a below the curve. A steeper 
than predicted drop in ISS with redshift could imply redshift evolution in the Doppler factor of the 
jets or angular broadening caused by propagation through the intergalactic medium, as mentioned 
by Lovell et al. (2008). However, the discrepancy shown here is not statistically strong enough 
and we await further redshift values from the remaining (weak sample) sources before further 
analysis. In fact, this discrepancy at the highest redshift bin could be linked to a steepening of the 
mean source spectral index with increasing redshift, coupled with the lower D(2days) of steeper 
spectrum sources, as has been found in the dual-frequency follow-up observations (Koay et al. 2012). 
As noted in Section 2, the flux densities used to estimate the spectral indices for the selection of 
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MASIV sources were non-simultaneous; and the dual-frequency follow-up observations provided 
more accurate estimates. Indeed any more detailed analysis must also consider the distributions in 
the parameters of the model rather than the single values adopted here. Nevertheless, we stress that 
the present work presents the strongest observational evidence to date for the predicted increase in 
angular size due to cosmological expansion. 


6. SUMMARY 

We present the optical and spectroscopic identifications of a sample of 475 flat-spectrum radio 
sources which are unresolved at 5 GHz (FWHM < 50 mas). About half of these sources also show 
ISS. Our sample is divided into radio strong (S $ghz > 0.3 Jy) and radio weak ( 0.06< S 5 ghz <0.3 
Jy) categories. For optical identification the SDSS DR5 and GSC2.3 surveys were used. About 
80% of the radio strong and 60% of the radio weak sources have an optical counterpart with R < 
20 magnitude. Using sub-samples defined from the SDSS data indicates that the nearly complete 
radio strong sample ((S^ghz) = 0.9 Jy, N=78) has median r-magnitude ~ 17.9 (r max ~ 22.0), and 
the 70%-identified radio weak sample has r ~ 18.9 (r max ~ 22.9, (Ssghz) = 0.1 Jy, N=95). 

Spectroscopic identifications and redshifts are from the literature and from our own observa- 
tions (347 sources). Spectroscopic identification of the radio strong sample is 91% complete and 
the weak sample is 57% complete. 

The key findings can be summarized as follows: 

1) The radio weak sources have flatter radio-optical spectral index than the radio strong 
sources. 

2) Almost 80% of the MASIV sources are identified as broad line objects, ~13% as BL Lacs 
and ~7% as narrow line objects or galaxies. The spectroscopic identifications are similar for the 
radio strong and weak samples. Our results suggest that selecting for compact radio structure 
favors FSRQs and BL Lacs over other flat spectrum sources. 

3) Of the scintillating sources, 25% are BL Lacs and the rest are almost exclusively broad 
line objects. The fraction of BL Lacs increases with the degree of persistence of ISS and 40% of 
the sources with most persistent ISS (scintillating at all four epochs) are BL Lac objects. This 
is consistent with BL Lacs being more strongly core dominated than FSRQs, and could possibly 
indicate relativistic jets that are pointed closer to the direction of the Earth. We confirm that 
D(2days) depends on 5 GHz flux density, with weak sources having higher D(2days) values. A 
similar trend is found also when broad line or weak line objects are studied separately. We found 
indications that the ISS properties correlate with emission line equivalent width. 

4) Radio and optical luminosities of the scintillating and non-scintillating broad line sources 
are similar at a given redshift, suggesting intrinsically similar SEDs. 
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5) We find no correlation between the optical luminosity and ISS. We find a weak correlation 
between radio power and ISS. However at given redshift the ISS and non-ISS sources have similar 
radio power. 

6) For broad line objects, we confirm the sharp decrease in the number of ISS sources and in 
the level of their ISS above a redshift ~ 2. The decrease is compared to a simple model for ISS 
of flat spectrum radio sources with maximum brightness temperatures that are Doppler-boosted 
in jets pointing toward the Earth. There is reasonable agreement but more redshifts are needed 
for the weak radio sample. We found strong observational evidence for the predicted increase in 
angular size due to cosmological expansion. 
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Fig. 2.- The histogram of the magnitude difference for objects with both SDSS and GSC identi- 
fication. The bottom panel shows the converted R-magnitude difference. 





Fig3-eps-converted-to .pdf 


Fig. 3. — 5 GHz flux density vs optical flux. The blue circles, red stars and open circles are 
for variable, non-variable and intermittently variable, respectively. The magnitude upper limits 
(R=20.5) are indicated with blue, red, black triangles for variable, non-variable and intermittently 
variable, respectively. The diagonal lines indicate constant radio-optical spectral index -0.2, -0.45, 
-0.55 and -1.0 from top to bottom. 





Fig. 4. — Example spectra of an FSRQ (Type 1), NLRG (Type 2) and BL Lac object. In the top 
panel the vertical bars indicate redshifted Lya, Si II, C IVA1549, C IIIA1909 and Mg IIA2798. In 
the middle panel vertical bars indicate redshifted Ne VA3426, 0 IIA3730, Ne IIIA3870, Ne IIIA3870, 
0 III]AA4958, 5007 and Ha. In the bottom panel the vertical bars indicate absorption lines at 
redshift 1.41 with Mg IA2852, Mg IIAA2803,2796. The ©-signs identify features due to the Earth’s 
atmosphere 
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Fig. 6. — Optical R-band luminosity versus radio power, with FSRQ, BL Lac and AGN marked 
with red circles, blue squares and green stars respectively. The filled symbols show in the top panel 
the location of optically bright sources, in the middle panel radio weak sources and in the bottom 
panel redshift 1.3 < z < 1.8 objects. Note the two outliers (J0248+0434 and J1141+5953) near 
Mr=- 26, whose magnitude estimates from the USNO-B1 are too bright. 
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Fig. 8. — D(2days) against redshift for the Type 1 sources. Red triangles indicate radio strong 
sources and blue circles radio weak sources. The horizontal thin and thick lines indicate the median 
Z)(2days) for the given redshift interval for radio strong and weak samples, respectively. Type 0 
objects are not shown as most of them do not have redshifts. 





Fig. 9.— The data for all FSRQs averaged into 4 redshift bins. Over plotted are theoretical curves 
for a simple ISS model, discussed in the text, with Doppler factors as indicated in the legend. The 
horizontal line is at 4 xlO -4 , which is the threshold for significant variation. 




Table 1. Summary of the optical identification using GSC2.3 and SDSS DR5. 


5 G tf 2 limit 

N 

median 

EF 

survey 

S>0.3Jy 

193 

18.2 

36 

GSC2.3 F 

S<0.3Jy 

167 

19.5 

79 

GSC2.3 F 

7 h 

< RA < 16 h and 5 < 64 

S>0.3Jy 

62 

18.25 

4 

SDSS r 

S<0.3Jy 

75 

19.82 

11 

SDSS r 


Note. — Columns:(l) the radio flux density 
limit, (2) the number of detected optical counter 
parts (3) the median magnitude of the sample (no 
galactic extinction correction) (4) non-detections 
(5) The survey and pass band 



A. COMMENTS ON ARCHIVAL REDSHIFTS 


After careful inspection three objects with archival redshift were rejected. In addition, un- 
reliable redshifts have been proposed for ten BL Lac objects and four objects with adopted, but 
uncertain redshift or spectroscopic identification. We briefly discuss these here. 

J0200+0322: Fricke et al. (1983) reported redshift, 2 = 0.765, based on a single (Mg IIA2798) line 
at 4920A. Their spectrum had the wavelength range from 4200-7100A where for a typical QSO 
spectrum a single broad line can be detected at z ~ 0.8 (Mg II) or 1.6 (C III]). 

J0204+1514: Stickel et al. (1996) reported 2=0.833 based on two narrow emission lines. Later 
Perlman et al. (1998) found more lines and a better fit for 2=0.405 and they classified this object 
as an FSRQ, which we adopted for this work. However this object should be re-observed in order 
to secure the spectroscopic identification. 

J0406+2511: Hook et al. (1996) assumed the single line at 6260A to be [O II] A3727, at 2=0.68. 
Their data has a wavelength range from 5200-9100A, however the red and blue ends are very noisy, 
reducing the useful range to 5800-8500A. If the redshift is indeed 0.68, H j3 should be around 8150A, 
and/or [O III]A5007 around 8410A, which are not seen. 

J0449+1121: von Montigny et al. (1995) reported 2=1.27, however Halpern et al. (2003) found 
featureless spectrum. Thus we adopted a BL Lac identification for this source. 

J0738+1742: NED reports 2 = 0.424. However, this is the redshift of intervening Mgll absorption 
feature, i.e. 2 > 0.43 (Rector & Stocke (2001) and references therein) 

J0745+1011: We adopted the redshift from Best et al. (2003) who measured 2=2.624, however, 
Labiano et al. (2007) could not confirm this. If the redshift is correct the radio power is 4x 10 28 W /Hz 
and an outlier in the MR-Log(.P)-plot (Figure 6). This object should be re-observed. 

J0753+5352: This source has featureless optical spectrum and the redshift has a lower limit of 0.2 
based on imaging data (Stickel & Kuehr 1993). 

J0818+4222: This is a radio source with featureless optical spectrum and has no reliable redshift 
available (Rector & Stocke (2001) and references therein). 

J0929+5013: The SDSS redshift is uncertain as there are no obvious emission lines. 

.11150+2417: This object has a featureless optical spectrum Rector & Stocke (2001, and references 
therein) . 

.11309+1154 This object has a featureless optical spectrum by Lynds & Wills (1972) and SDSS. 
Note that this object has SDSS DR5 redshift of 0.3852 and SDSS DR 8 redshift of 2.601. 
.11446+1721: The redshift of this FSRQ is uncertain as Sowards-Emmerd et al. (2005) reported, 
2=0.102, however Healey et al. (2008) reported 2=1.026. For this work we adopt FSRQ classifica- 
tion with unknown redshift. 

J1502+3350: White et al. (2000) report featureless spectrum for this object. Note that the SDSS 
DR6 suggest 2 = 2.222, however the spectrum appears featureless. 

J1602+3326: The redshift is based on optical magnitude (Snellen et al. 2000) using the GPS galaxy 
Hubble diagram (Stanghellini et al. 2005) 

J1648+2141 The SDSS spectrum is featureless, however the SDSS DR6 suggest 2=1.085. 



Table 2. MASIV Data 
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Table 3. Objects with single emission line from our data. The range indicates the useful 

wavelength range of the spectrogram. 


Source 

Line 

A 

Range 

A 

S/N max 
per pixel 

J0213+3652 

4179 

3500-5000 

2 

J0342+3859 

5453 

4150-7500 

15 

J1316+6927 

5675 

4500-7500 

3 

J1711+6853 

5660 

3900-8500 

8 

J1953+3537 

5125 

5200-7500 

10 

J2012+5308 

3967 

3950-7500 

6 


Table 4. Summary of the spectroscopic identification. The third and fourth column indicates the 
number of sources in the radio strong and radio weak samples, respectively. 


Type 

Spect id 

Strong 

Weak 

Type 0 

BLL 

37 

11 


BLL? 


5 


PEG 


6 

Type 1 

FSRQ 

155 

107 


FSRQ? 

1 



Sy 1 


2 


Sy 1/BLRG 

1 



NLSyl 

1 



NLQSO 

2 



Syl. 5/LINER 

1 


Type 2 

Sy2 

3 

2 


Sy2? 


1 


NLRG 

5 

1 


LINER 

1 

1 


HU? 


1 


galaxy 


3 

Single line objects 

1 

5 
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J1649+7442: This redshift is from Appenzeller et al. (1998) from the spectrum with a wavelength 
range from 4000 to 7000/9000A. The redshift is based on a single line at 6658A which is assumed 
to be Mg IIA2798 at 2=1.378. At this redshift, one would expect to see C III] A1909 at 4532A, 
however if the line was C III] or H f5 one would expect to see C IV at rsj 5392A or [O III] near 6870A. 
The redshift is adopted, but needs to be confirmed. 

J1719+1745: Sowards-Emmerd et al. (2005), reported 2=0.173. However Shaw et al. (2009) find 
featureless spectrum and they estimate redshift to be >0.58. This source is classified as a BL Lac 
object without known redshift. 


B. TRANSFORMATION FROM GSC2.3 TO LANDOLT SYSTEM 

In order to transform the GSF2.3 magnitudes to the Landolt (2009) system, we obtained 
JFN magnitudes from the Landolt standard stars. We included only stars with four or more 
measurements by Landolt, then plotted J — B, F — R and N — I and removed outliers. Least 
squares polynomial fit resulted following transformations: 


Base = J + 0.18 (J - F) — 0.08, 

(Bl) 

Rose = B[~ 0.085(.F — N ) + 0.04], 

(B2) 

Iqsc = 1V[+0.066(F — N) — 0.02], 

(B3) 


The terms in square brackets for Rose and lose are from the best least-squares fit, however 
the fit is improved only marginally, hence we have omitted the color term. The root mean squared 
error for the BRI-conversions are 0.10, 0.11 and 0.11, respectively. 

The Base transformation is valid for the ( J — F)-colors from -0.5 to 2.5 and Rose, lose for 
the (F — A)-color range from -0.4 to 1.4. 

C. NOTES ON INDIVIDUAL OBJECTS 

In this appendix, we discuss some of our spectroscopic follow-up targets. For the narrow line 
objects we report the EW of the strongest line. 

J0019+2021 BL Lac classification is from Chu et al. (1986) who found no strong emission lines. 
Our DBSP data (S/N<35) shows absorption system at 2=1.41 with Mg IA2852, Mg IIAA2803,2796, 
(Fe IIA2344) and a possible absorption system at 2 = 2.274 (Figure 4). 
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J0150+2646 The DBSP (07-2007) and NOT (8-2004) spectrum are featureless with S/N<14 and 
S/N<8, respectively. Thus, this source is tentatively identified as a BL Lac object. 

J0509+0541 Halpern et al. (2003) reported featureless spectrum from several observations. Our 
high S/N (>130) data is featureless showing only galactic absorption lines and unidentified absorp- 
tion at 4276A 

J0800+4854 Our data show strong absorption at the C IVA1549- and Lya-lines and also the 
He IlA1640-line is easily detected. 

J1 118+2922 This low redshift object has two easily detected narrow emission lines at A3998A and 
A7038A, identified as [O II] A3727 and Ha at redshift of 0.072. The line ratios suggest HII region 
or Liner classification. 

J1352+3603 We detect narrow [O III] AA4958, 5007 and [O II] A3727 lines at z=0.3057 with, M B =- 
21.5, (Mr=- 23.0), C=42%, EW([0 II])=20A and noisy EW(Ha)~ 37A. The source is classified as 
PEG based on the Ca-break and line EW. On the other hand, based on the classification scheme 
of Landt et al. (2004), this is a strong-lined source. 

J1442+0625 Our data show a single line at 4753 A, with continuum from 4800-7500A (S/N<10). 
Our tentative identification is Mg IlA2798-line at z=0.698. If this line would be C III] at z=1.49 
we would expect to see Mg IIA2798 at about 7000A which is not detected. 

J1444+0257 This object has a featureless spectrum from 4500 to 7500A with S/N<5, possibly a 
BL Lac object 

J1444+1632 We detect a single line at 5688A with noisy continuum (S/N <3) and tentative line 
at 7247A, these are identified as Lya and C IVA1549, which would indicate z=3.7. 

.11505+0326 This object has been classified as FSRQ in CGRaBS (Healey et al. 2008). Our NOT 
data confirms the CGRaBS and SDSS redshift. However reanalyzing the SDSS data and based on 
the NOT data the rest frame line widths are less than 2000km/s indicating NLSyl classification. 
.J1651+0129 Our DBSP data shows a possible line at 7128A with the continuum S/N<7. As we 
are not able to identify the line or spectrum, this target is not included to the analysis. 
J1718+4448 The DBSP red spectrum show no obvious lines from 5300 to 8300A with S/N < 5, 
possibly a BL Lac object. 

J1728+1931 We detect narrow ([O II] A3727) , [O III] AA4958, 5007 and a noisy [O I]A6300 lines at 
z=0.1759, with [O III]A5007/H/3>3, EW([0 II])=13A, M B =-21.1, (M R =-24.2) and C=42%. The 
Ha-line is noisy, however it appears to be narrow. This is border line object between strong and 
weak lined sources using Landt et al. (2004) classification. We tentatively classify this as a Sy2-type 
object. 

J1734+3857 Stickel et al. (1989) classification for this source is an FSRQ at z=0.97 based on a 
single broad emission line which was assumed to be Mg IIA2798 (EW(Mg II)=15.3A, observed) Our 
data shows also only a single line at 5535A assumed to be Mg II at ,2=0.97 with the EW=4.2A, 
hence BL Lac classification. 

J1 739+4737 This is a BL Lac-object with featureless spectrum (Scarpa et al. 2000). Our DBSP 
data show no lines with S/N<20 

.11742+5945 This object has a featureless spectrum Laurent-Muehleisen et al. (1998), and the red- 
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shift (z=OA) is based on a detection of the host galaxy (Nilsson et al. 2003). Our NOT spectrum 
(07-2006) shows no lines from 4000 to 7500 A with S/N<25. 

J1 747+4658 This source has no emission line redshift. Vermeulen et al. (1996) reported an inter- 
vening Mg II absorption line system at z= 1.484, which, by some authors, has been adopted as the 
source redshift. Our NOT spectrum (Grisrn 6) from 4000 to 5500A with S/N<6 is featureless. 
J1757+0531 Narrow emission lines, [O II]A3727 and [O III] AA4958, 5007 at z= 0.34482 are clearly 
visible and so is a noisy strong Ha feature at the red end of the spectrogram. The line and continuum 
properties suggest Type 2 identification: [O IIl]A5007/H/3>3, Mr=- 23.2, EW([0 III]A5007)=37A and 
C=28%. However the continuum is rather noisy, (S/N<8) and weak broad emission lines might 
have been undetected (Figure 4). Tentative spectroscopic identification is NLRG, however, based 
on Fig. 8 and 9 of Caccianiga et al. (2008) this object could be a Type 1 object with the AGN core 
diluted by the host galaxy. 

J1832+1357 Our data confirms the CGRaBS (Healey et al. 2008) redshift. We detect also some 
absorption features in C IVA 1549- line 

J1848+3219 Sowards-Emmerd et al. (2005) reported this source as an FSRQ at 2=0.798. Our 
data shows a strong line at 5035A, a noisy narrow line at 6703A (identified as Mg IIA2798 and 
[O II] A3727, respectively) and no strong lines in 3600-5000 and 5100-7500A with S/N<17, sup- 
porting the literature identification and redshift. 

J1905+1943 A strong broad line is detected at 5106A and a tentative broad line at 6267A, sug- 
gesting z =2.3. 

J1950+0807 From the DBSP data we detect narrow [O II] A3727, [O III]AA4958, 5007 and broad 
Ha-lines at 2=0.2964, with Mr=-22.8, Mr=- 20.7 and EW([0 III]A5007)=-48A, suggesting Syl.9 
classification. 

J2023+5427 The redshift 1.48 is measured from Mg IIA2798 and C III] A 1909-lines. In addition 
we detect Mgl(2852), MgII(2803,2796), FeII(2387,2600) and FeII(2383,2344) absorption lines at 
2=1.415 

J2130+0339 From the DBSP data we detect narrow lines [O II] A3727 and [O III] AA4958, 5007 at 
2=0.653, with EW([0 III]A5007)=140A, Mr >-22.5 and C=31%, suggesting Type 2 identification. 
Note, this source is not detected in the GSC2.3. 

J2155+0916 This object is tentatively identified as a BL Lac object. Our NOT spectrum is fea- 
tureless from 4500 to 7500A with S/N<20, however if the redshift is ~0.6, both Mg IIA2798 and 
H/I-lines will be outside the spectrogram. 

J2201+5048 We confirm the Sowards-Emmerd et al. (2005) redshift, 1.899 and in addition we detect 
Mg II-absorption system at 2=1.3822 with Mg IA2852, Mg IIAA2803,2796, Fe IIA2344, 2374, 2383, 
however this source was rejected from the Lovell et al. (2008) analysis due to evidence of source 
structure, hence not included to the analysis. 

J2203+1725 Smith et al. (1977) classified this source as a FSRQ based on strong C III]A1909- and 
Mg IlA2798-lines at 2=1.076. Our data shows a single line at 581lA with EW(Mg II)=3.3A and 
S/N<60 around the line, hence BL Lac object classification. 

J2212+2759 This is a BL Lac object, with no detected emission lines, based on our DBSP data 



-45 - 


(07-2007) with S/N<45 and NOT data (07-2005) with S/N<18. 

J2230+6946 We confirm the Healey et al. (2008) (CGRaBS) redshift. 

J2237+4216 An emission line is detected at 5655A and this is the only emission line between 4500 
and 7500A, suggesting Mg II identification. 

J2241+4120 Henstock et al. (1997) and Sowards-Emmerd et al. (2005) classified this as a BL Lac 
object. Our DBSP data (07-2007) with S/N<60 and NOT data (07-2005) from 4000 to 7500A with 
S/N<30 confirm this. 

J2242+2955 A strong broad line is detected at 4102A and a tentative broad line at 5013A, sug- 
gesting z =1.7. 

J2258+0516 We detected a strong absorption feature at the centre of the C IVA1549-line and Lya 
line is almost fully absorbed. 

J2303+1431 This is a BL Lac object with featureless spectrum. Our DBSP spectrum (07-2007) 
has S/N<25 and NOT spectrum (07-2006) S/N<8. 

J2315+8631 A broad line at 5273A (Mg IIA2798) is visible in our DBSP (08-2007) and NOT 
(07-2006) spectrum. This is the only strong emission line between 3700 and 8500A (S/N<30) , 
hence a tentative redshift of z=0.88. 

J2325+3957 Shaw et al. (2009) reported a featureless spectrum with intervening Mg IIA2798 and 
Fell absorption features at z= 1.04. Our DBSP data (07-2007) shows no emission lines. However, 
we detect absorption features which are identified as Mg IA2852, Mg IIAA2803,2796, Fe IIA2600 
(and Fe IIA2344) at z=0.4695 and Mg IIAA2803,2796, Fe IIA2600 at z=0.4156 with the continuum 
S/N up to ~25. We also detected an absorption feature at 5757Apossibly Mg IIA2798 at 2 = 1.04, 
but not the expected Fe II feature at 5353A. The 4119A feature, Mg IIA2803 at z=0.4695, is also 
visible from the NOT spectrum (07-2006). 
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Table 5. Classification break down of the radio strong and weak samples including 
calibrators. The ’SDSS sample’ source distributions are in column four and 


the secondary 
five. 


type N full % N S dss % 


Radio strong sources 


Type 0 

37 

184=3% 

19 

244=6 % 

Type 1 

162 

784=6% 

57 

724=10% 

Type 2 

9 

44=1% 

2 

34=2 % 

Radio weak sample 

Type 0 

22 

164=3% 

8 

144=5% 

Type 1 

108 

784=7% 

46 

784=11% 

Type 2 

9 

64=2% 

5 

84=4% 


Note. — The errors are Poissonian. 



-47« 


Table 6. Summary of some flat radio spectrum AGN surveys. Note, the PEGs in CBS and 
MASIV are included in the “other AGN” as is the case with CGRaBS and DXRBS. 


Survey 

5 GHz limit 
mJy 

R-rnag 

FSRQ 

% 

BL Lac 

% 

other AGN 

% 

Comment 

CBS 

30 

17.5 

45 

21 

34 


CBS rb 

65: 

17.5 

51 

20 

29 

MASIV 5 GHz limit 

CGRaBS 

65 

23.7 

88 

10 

2 

Selection by Figure of M erit 

DXRBS 

50 

24.0 

80 

12 

8 

some steep spectrum sources 

MASIV 

60 

23.1 

78 

15 

7 

Selection by 8.4 GHz 

MASIV ob 

60 

17.5: 

64 

24 

12 

Optically bright CBS limit 

MASIV ISS 

60 

21.3 

64 

34 

2 

Persistent IDV sources, IDV3, 4 

MASIV SDSS 

60 

21.7 

75 

18 

7 

SDSS sample 


Note. — Columns: (1) the name of the sample, (2) 5 GHz radio flux density limit for faintest objects, (3) 
optical magnitude limit or the faintest detected objects, (4) the fraction of broad line objects, (5) the fraction 
of BL Lac objects (6) the fraction of other spectral type objects, mostly narrow line objects. 

The colon indicates MASIV radio flux density limit for the “CBS rb” , and CBS optical limit for “MASIV ob” . 
Note, the DXRBS includes steep radio spectrum sources as well (Landt et al. 2001). 


Table 7. Median D(2days) of the subsamples selected by optical spectral type and 5 GHz flux 
density. The First column indicates optical spectral type, followed by the median D(2days) and 

the number of objects. 


Type 

D(2days)xl000 

N 

D(2days)xl000 

N 


S^ghz > 0.3 Jy 


Ssghz < 0.3 Jy 


0 

0.701 

31 

2.35 

22 

1 

0.335 

143 

0.868 

108 

2 

0.193 

6 

0.44 

9 


Type 1 z 

< 1.4 

sources 


1 

0.346 

79 

1.27 

54 
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Table 8. Number of sources and the fractions of “variable” MASIV epochs for the three optical 
spectral types. The first column indicates the number variable in epochs and the next columns 
show the fraction of sources in each group with the number of sources in parentheses. 


Sample 

IDV 

Type 0 

strong weak 

Typ' 

strong 

e 1 
weak 

Type 2 
strong weak 

0 

13% (4) 

18% (4) 

36% (51) 

36% (39) 

5 7 

1 

19% (6) 

18% (4) 

24% (34) 

15%(16) 

1 1 

2 

13% (4) 

(0) 

24% (34) 

17%(18) 

0 0 

3 

29% (9) 

18% (4) 

8% (12) 

18%(19) 

0 1 

4 

26% (8) 

46% (10) 

8% (12) 

15%(16) 

0 1 



Table 9. Redshift, core radio power and optical luminosity, excluding secondary calibrators. The 
IDVO, IDV1, etc. indicates the number of epochs of variability (Lovell et al. 2008). 


Strong sources Weak sources 


Sample 

N 

< z > 

log (P r ) 

Wm~ 2 

Mr 

N 

< z > 

log(-P r ) 

Wm~ 2 

Mr 

Type 0 

16 

0.4T0.3 

26.6=40.6 

-25.2=41.1 

10 

0.2=40.1 

25.0=40.8 

-23.6=41.3 

Type 1 

142 

1.3T0.7 

27.6=40.6 

-26.4=41.6 

108 

1.4=40.8 

26.8=40.5 

-25.9=41.6 

Type 2 

6 

0.4T0.3 

26.7=41.3 


9 

0.1=40.2 

24.1=41.1 






Type 1 





.D(2days)> 4 x 10~ 4 

59 

1.24=0.6 

27.5=40.6 

-26.2=41.4 

76 

1.3=40.7 

26.7=40.5 

-25.7=41.5 

D(2days)< 4 x 10" 4 

83 

1. 4=1=0. 7 

27.7=40.5 

-26.7=41.6 

32 

1.8=41.0 

27.0=40.5 

-26.7=41.4 





1.0 < z < 1.8 





.D(2days)> 4 x 10 -4 

24 

1. 4=1=0. 2 

27.5=40.2 

-26.4=41.0 

29 

1.3=40.2 

26.7=40.2 

-26.0=41.2 

.D(2days)< 4 x 10 -4 

35 

1. 4=1=0. 2 

27.6=40.2 

-26.4=41.3 

10 

1.3=40.3 

26.8=40.2 

-25.7=41.0 





Type 1 





ISS (IDV234) 

58 

1. 2=1=0. 7 

27.5=40.6 

-26.3=41.4 

53 

1.2=40.6 

26.7=40.5 

-25.4=41.4 

IDVO 

50 

1. 4=1=0. 8 

27.7=40.6 

-26.7=41.8 

39 

1.9=41.0 

27.1=40.5 

-26.5=41.6 

IDV1 

34 

1.3=40.6 

27.5=40.5 

-26.3=41.3 

16 

1.4=40.7 

26.8=40.4 

-26.0=41.3 

IDV2 

34 

1.5=40.7 

27.6=40.6 

-26.4=41.5 

18 

1.3=40.7 

26.7=40.5 

-26. Oil. 3 

IDV3 

12 

1.0=40.5 

27.4=40.6 

-25.9=41.2 

19 

1.1=40.7 

26.6=40.6 

-25. Oil. 4 

IDV4 

12 

1.4=40.7 

27.6=40.4 

-26.5=41.3 

16 

1.2=40.5 

26.7=40.4 

-25.3il.4 


Note. — Columns: (1) the sub-sample, (2) Number of objects, (3) Median redshift, (4) Median of 
the Log radio power with standard deviation, (5) Median absolute magnitude with standard deviation, 
(6)-(9) same as Cols. (2)-(5). 



